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^ ; Abstract 

This paper investigates the physics reach of the IceCube neutrino detector when it will have 
collected a data set of order one million atmospheric neutrinos with energies in the 0.1 ~ 10 4 TeV 



range. The paper consists of three parts. We first demonstrate how to simulate the detector 
performance using relatively simple analytic methods. Because of the high energies of the neutrinos, 
their oscillations, propagation in the Earth and regeneration due to r decay must be treated in a 
coherent way. We set up the formalism to do this and discuss the implications. In a final section 
we apply the methods developed to evaluate the potential of IceCube to study new physics beyond 
neutrino oscillations. Not surprisingly, because of the increased energy and statistics over present 
experiments, existing bounds on violations of the equivalence principle and of Lorentz invariance 
can be improved by over two orders of magnitude. The methods developed can be readily applied 
to other non-conventional physics associated with neutrinos. 
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I. INTRODUCTION 



After the discovery of neutrino oscillations in underground experiments, the observations 
have been confirmed by experiments using "man-made" neutrinos from accelerators and 
nuclear reactors 0|. We are entering an era of precision neutrino physics. In this context 
we discuss the unique potential of IceCube, an experiment that will collect large statistics 
samples of high energy atmospheric neutrinos. In contrast with its other missions, the beam 
and its physics exploitation are guaranteed. 

With its high statistics data [2j Super-Kamiokande (SK) established beyond doubt that 
the observed deficit in the /z-like atmospheric events is due to oscillations, a result also 
supported by the KEK to Kamioka long-baseline neutrino oscillation experiment (K2K) 0] 
and by the MACRO and Soudan 2 p experiments. 

It has been recognized that oscillations are not the only possible mechanism for atmo- 
spheric — > v T flavour transitions [7] . These can also be generated by a variety of nonstan- 
dard neutrino interactions characterized by the presence of an unconventional interaction 
(other than the neutrino mass terms) that mixes neutrino flavours [jj. Examples include 
violations of the equivalence principle (VEP) @,0,EI3|> non-standard neutrino interactions 
with matter nj], neutrino couplings to space-time torsion fields 3], violations of Lorentz 
invariance (VLI) (lH Q] and of CPT symmetry (15I 16 . 17]. From the point of view of 
neutrino oscillation phenomenology, a critical feature of these scenarios is a departure from 
the E energy dependence of the conventional oscillation wavelength 18 , 3] . Although these 
scenarios no longer explain the data|2ol. 21, 22, 23, 24, 25], a combined analysis of the at- 



mospheric neutrino and K2K data can be per 
on the size of subdominant oscillation effects 



brmed to obtain the best constraints to date 
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In contrast to the E energy dependence of the conventional oscillation length, new physics 
predicts neutrino oscillations with wavelengths that are constant or decrease with energy. 
IceCube, with energy reach in the 0.1 ~ 10 4 TeV range for atmospheric neutrinos, is the 
ideal experiment to search for new physics. For most of this energy interval standard Am 2 
oscillations are suppressed and therefore the observation of an angular distortion of the 
atmospheric neutrino flux or its energy dependence provide a clear signature for the presence 
of new physics mixing neutrino flavours. 

In this paper we explore the physics that can be probed with the high-statistics high- 
energy atmospheric data that will be collected by the IceCube detector. In particular we 
quantify its sensitivity to atmospheric neutrino oscillations driven by new physics effects. 
The outline is as follows: Our analytic "simulation" of the IceCube detector is described 
in Sec. |H] where the expected number of atmospheric neutrino events and their energy 



1 Atmospheric neutrinos have also been used to place bounds on other exotic forms of new physics in the 
neutrino sector such as the possibility of neutrino decay |27l bgf or quantum decoherence in the neutrino 
ensemble El El 
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distribution are presented. In Sec. II I II we briefly summarize the formalism for discussing the 
phenomenolgy of non-standard neutrino oscillations and we derive the evolution equations 
that describe a high-energy neutrino beam subject to oscillations as well as attenuation and 
v regeneration due to r decay In Sec. IIVI we illustrate the sensitivity of the 

detector for violations of Lorentz invariance and the equivalence principle. 



II. SIMULATION OF MUON EVENT RATES IN ICECUBE 

In a high energy neutrino telescope muon neutrinos are detected via their charged current 
(CC) interactions in the matter surrounding the detector. Such interactions produce muons 
which reach the detector. High energy muons have very large average range resulting in an 
effective volume of the detector significantly larger than the instrumented volume. 

In our semianalytical calculation we will obtain the expected number of induced events 
from 

/l /*oo /»oo /*oo poo 

dcos6 dl dEf / dEl / <LE V (1) 
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,,. , a is the differential muon neutrino neutrino flux in the vicinity of the detector after 

dE v a cos J 

evolution in the Earth matter (see next section for details). We use as input the neutrino 



fluxes from Honda [35] which we extrapolate to match at higher energies the fluxes from 
Volkova 36]. At high energy prompt neutrinos from charm decay are important. In order to 
estimate the uncertainty associated with the poorly known charm meson production cross 
sections at the relevant energies, we compute the expected number of events for two different 
models of charm production: the recombination quark parton model (RQPM) developed by 
Bugaev et al (3?| and the model of Thunman et al (TIG) ji^ that predicts a smaller rate. 
-^&(E U , E®) is the differential CC interaction cross section producing a muon of energy E° 
and tit is the number density of nucleons in the matter surrounding the detector and T is 
the exposure time of the detector. Equivalently, muon events arise from interactions that 
are evaluated by an equation similar to Eq.(J2J). 

After production with energy E®, the muon ranges out in the rock and in the ice sur- 
rounding the detector and looses energy. We denote by F(E^Ef t n ,l) the function that 
describes the energy spectrum of the muons arriving at the detector. Thus F(E%, E^ 1 , 1) 
represents the probability that a muon produced with energy E® arrives at the detector 
with energy E^ n after traveling a distance /. We compute the function F(E®, E^ n , I) by 
propagating the muons to the detector taking into account energy losses due to ionization, 
bremsstrahlung, e + e~ pair production and nnc.ear interactions according to Ref. 0. We 
include in F(E®, E^ n , I) the possibility of fluctuations around the average muon energy loss 
(using the average energy loss would equalize I to the average muon range distance). Thus 



3 



in our calculation we keep E®, E^ n , and, / as independent variables. For simplicity we use 
ny and F(E®, E^ n , I) in ice and we account for the effect of the rock bed below the ice in 
the form of an additional angular dependence of the effective area for upward going events 
(see Eq. JTUl below)). 

The details of the detector are encoded in the effective area A\^. We use the following 
phenomenological parametrization of the A Q e ^ to simulate the response of the IceCube de- 
tector after events that are not neutrinos have been rejected (this is achieved by quality cuts 
referred to as "level 2" cuts in Ref. [loT ]) 

A° eff = A (E f ; n ) x R(cos6, E^) x R(l mm ) . (2) 

In Aoi^Ejj 1 " 1 ) we include the energy dependence of the effective area due to trigger require- 
ments (see Ref. ^(J). We find good agreement with the results for the experiment MC 
simulation if we introduce a simple straight line dependence on log 1Q (E^ 



A (Ep n ) = A 



1 + 0.55 log 



10 



GeV 



(3) 



where Aq is an overall normalization constant which is fixed to reproduce the expected num- 



ber of events in the absence of oscillations: 91000 events/yr after level 2 cuts 



'or conventional 



atmospheric neutrinos. We next have to "simulate" cuts introduced in Ref. 4(J in the niuon 
tracklength l min and the number of optical modules reporting signals in an event NcH,min- 
Ri}min) represents the smearing in the minimum track length cut, l min = 300 m, due to 
the uncertainty in the track length reconstruction which we parametrize by a Gaussian 



(7 • -V 

\<"min l m 



\2 



RiJ-min) — r- — / dl min exp „ , (4) 



2nai Jo 2a 



with G\ = 50 m. 

The angular dependence of the effective area for downgoing events (9 < 80°) is determined 
by the level 2 cut on the minimum number of channels Nqh > NcH,min{ c °s@) — 150 + 
250 cos#. We translate this requirement in an E^ m -dependent angular constraint as 

R(<x*0,E}r) = -= / dN CH exp — 2 ^— , (5) 

,min 

(6) 

where (Nch) E fin ^ s ^ e average channel multiplicity produced by a muon which reaches 
the detector with energy E^ n and cr^ Cfr is the spread on the distribution. Using Fig. 7 in 
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Ref. 40] we obtain the parametrization 

log 10 ((iVaff)^) = 2.0 + 0.88 

X = 0.47|log 10 M-4.6 



X 



a NcH = 0.4 (N CH ) 



(7) 

(8) 
(9) 



Finally, we can account for the presence of the rock bed below the detector by introducing 
a phenomenological angular dependence of the effective area for upward going muons 



i?(cos 9) = 0.70 - 0.48 cos 9 for 9 > 85° 



(10) 



independent of the muon energy. 

We show in Fig. [I] the effective area A e ff, defined as the ratio of the number of upgoing 
muon events, with/without the inclusion of A (E^ n ) x R(l min ) and the level 2 cuts on i min , 
and compare our results to the detector simulations after cuts from Fig. 5 of Ref. 0|. Our 
calculation correctly reproduces the energy threshold of the effective area. 
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FIG. 1: Effective function of the final muon energy after level 2 cuts in our calculation 

(full line) compared to the experimental MC simulation (data points). For comparison we also 
show ylo(-E> m ) (dashed line). 

Fig. |21 compares the energy spectrum and the zenith angular distribution of the events in 
the absence of oscillations after level 2 cuts obtained from our calculation with the results of 
the experimental MC. In both cases prompt neutrinos are included according to the RQPM 
model for charm production. The figure illustrates how our simple semianalytical calculation 
correctly reproduces the experimental simulation. 
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FIG. 2: Spectrum and zenith angular distribution after level 2 cuts for one year exposure obtained 
from our calculation (full lines) and from the experimental MC (dashed lines) (taken from Figs. 2 
and 9 of Ref. Q)- 

In Fig.Elwe show the expected spectrum of events in the absence of oscillations after level 
2 cuts as a function of the muon energy at the detector, E^ m (full line) . For comparison we 
also show the spectrum as a function of the muon energy before ranging E°. From the figure 
we read that in 10 years of operation IceCube will collect more than 7 x 10 5 atmospheric 
neutrino events with energies Ef in > 100 GeV. These events are generated by neutrinos with 
large enough energy for the standard Am 2 oscillations to be very much suppressed so they 
should behave as flavour eigenstates. This high-statistics high-energy event sample offers a 
unique opportunity to test new physics mechanisms for leptonic flavour mixing as we discuss 
next. 



III. PROPAGATION IN MATTER OF HIGH ENERGY 
OSCILLATING NEUTRINOS 

As described in the introduction, new physics (NP) scenarios can result in lepton flavour 
mixing in addition to "standard" Am 2 oscillations. We concentrate on v^—Vr flavour mixing 
mechanisms for which the pro pag ation of neutrinos (+) and antineutrinos (— ) is governed 
by the following Hamiltonian [161 ] : 

i A<5 n E n f-1 o\ TTt 

where Am 2 is the mass-squared difference between the two neutrino mass eigenstates, 
accounts for a possible relative sign of the NP effects between neutrinos and antineutrinos 
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FIG. 3: Expected spectrum of events in the absence of oscillations after level 2 cuts for one year 
exposure as a function of the muon energy at the detector, (full line). For comparison we also 
show the spectrum as a function of the initial muon energy E®. 

and A5 n parametrizes the size of the NP terms. The matrices XJe and ^J^ n) ± Vn are given by: 

/ cos# smd\ ( cos£ n sin£ n e ±ir ?"\ 

\ — sm 9 cos 91 V - sin £„e Tir > n cos £ n / 

by ?7„ we denote the possible non-vanishing relative phases. 

If NP strength is constant along the neutrino trajectory the oscillation probabilities take 
the form jlfi| : 

P v ^u, = 1 - = 1 - sin 2 26 sin 2 ( 1l) , (13) 



with 



sin 2 26 = — ^ (sin 2 26 + R 2 n sin 2 2£„ + 2R n sin 26 sin 2£ n cos r? n ) , ( 14) 



K = v 7 ! + Rl + 2R n (cos 26 cos 2£ n + sin 26 sin 2£„ cos 7] n ) , (15) 
,A5 n E n AE , . 

fl " = ff »^A^' < 16) 

where, for simplicity, we have assumed scenarios with one NP source characterized by a 
unique A5 n . 
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Eq. (fTTj) describes, for example, flavour mixing due to new tensor-like interactions for 
which n — 1 leading to a contribution to the oscillation wavelength inversely proportional to 
the neutrino energy. This is the case for u^s and z/ r 's of different masses in the presence of 
violation of the equivalence principle due to non- universal coupling of the neutrinos, 71 7^ 72 
(Vi and y?. being related to and v T by a rotation £„ ep ), to the local gravitational potential 

O0.0-EK 

For constant potential 0, this mechanism is phenomenologically equivalent to the break- 
down of Lorentz invariance resulting from different asymptotic values of the velocity of the 
neutrinos, C\ 7^ c 2 , with z/ x and u 2 being related to and v T by a rotation £ v n [13j,ll4 . 

For vector-like interactions, n = 0, the oscillation wavelength is energy-independent. This 
may arise, for instance, from a non-universal coupling of the neutrinos, k\ 7^ k 2 {v\ and v 2 
is related to the and v T by a rotation £q), to a space-time torsion field Q [l2j. Violation 
of CPT resulting from Lorentz- violating effects such as the operator, pfpfP^nVi -. a l so leads 



to an energy independent contribution to the oscillation wavelength 15, |16|, |13 which is a 
function of the eigenvalues of the Lorentz violating CPT-odd operator, hi, and the rotation 
angle, ^pt, between the corresponding neutrino eigenstates ^ and the flavour eigenstates 
v a . 

The flavour oscillations of atmospheric v^s in this scenarios is described by Eq. (fT3*|) with 
the identification: 

^ = Up A(5 1 = 2|0|(7 1 -7 2 )=2|0|A 7 <1.6xlO- 24 , for VEP (17) 
£i = &k, A6t = (ci -c 2 ) = 6c/c< 1.6 xlO" 24 , for VLI (18) 
^0 = , A5 = Q(ifei - k 2 ) < 6.3 x 10" 23 GeV , for coupling to torsion (19) 
£ = f^pT , A5 = h-b 2 < 5.0 x 10" 23 GeV , for 0PT , VLI (20) 

where for the first three scenarios a + = a~ while for the CPT violating case <r + = — a~ . 

At present the strongest limits on NP neutrino oscillations arise from the non-observation 
of departure from the Am 2 oscillation behaviour in atmospheric neutrinos at SK and the 
confirmation of oscillations with the same oscillation parameters from K2K. In Eqs. ()17f - 
1201) we quote the 3<r bounds from the up-to-date combined analysis of SK and K2K data 
performed in Ref. [2^. 

For most of the neutrino energies considered here, the standard Am 2 oscillations are 
suppressed and the NP effect is directly observed. As a consequence, the results will be 
independent of the phase i] n and we can chose the NP parameters in the range 

AS n >0, 0<£ n <7r/4. (21) 

The Hamiltonian of Eq. (jllj) describes the coherent evolution of the v^—Vr ensemble for 
any neutrino energy. High-energy neutrinos propagating in the Earth can also interact 



2 VEP for massive neutrinos due to quantum effects discussed in Ref. [Io| can also be parametrized as 
Eq. (HTJ with n = 2. 
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inelastically with the Earth matter either by charged current and neutral current (NC) and 
as a consequence the neutrino flux is attenuated. This attenuation is qualitatively and 
quantitatively different for i/ T 's and z/^'s. Muon neutrinos are absorbed by CC interactions 
while tau neutrinos are regenerated because they produce a r that decays into another tau 
neutrino before losing energy [3^. As a consequence, for each v T lost in CC interactions, 
another v T appears (degraded in energy) from the r decay and the Earth never becomes 
opaque to v' T s. Furthermore, as pointed out in Ref. (3^, a new secondary flux of z/^'s is also 
generated in the leptonic decay r — > pv^Vr. 

Attenuation and regeneration effects of incoherent neutrino fluxes can be consistently de- 
scribed by a set of coupled partial integro-differential cascade equations (see for example js^J 
and references therein). In this way, for example, the observed and oscillation-induced 
v T fluxes (and the associated event rates in a high energy neutrino telescope) from astro- 
physical sources has been evaluated. Alternatively, these effects can be accounted for in 



33, 



41 



Whatever 



a Monte Carlo simulation of the neutrino propagation in matter |3J 
the technique used, because of the long distance traveled by the neutrinos from the source, 
the oscillations average out and the neutrinos arriving at the Earth can be treated as an 
incoherent superposition of mass eigenstates. 

For atmospheric neutrinos this is not the case because oscillation, attenuation, and re- 
generation effects occur simultaneously when the neutrino beam travels across the Earth's 
matter. For the phenomenological analysis of conventional neutrino oscillations this fact can 
be ignored because the neutrino energies covered by current experiments are low enough for 
attenuation and regeneration effects to be negligible. Especially for non-standard scenario 
oscillations, future experiments probe high-energy neutrinos for which the attenuation and 
regeneration effects have to be accounted for simultaneously. In order to do so, we modify 
the neutrino flavour oscillation equations and couple them to the r evolution equations as 
described next. 

We find it convenient to use the density matrix formalism to describe neutrino flavour 
oscillations. The evolution of the neutrino ensemble is determined by the Liouville equation 
for the density matrix p(t) = u{t) ® u{t)^ 

f = "*M, (22) 

where H is given by Eq. (fTTj) . The survival probability in Eq. (|T3*|) is given by -P^(t) = 
Tr [n„ p{t)}, where U Ufi = ® is the z/ M state projector, and with initial condition 
p(0) = EL. An equivalent equation can be written for the antineutrino density matrix. 

For the case of oscillations between two neutrino states the hermitian operators p, H and 
the flavour projectors EL and EL r can be expanded in the basis formed by the unit matrix 
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and the three Pauli matrices 0{. In particular we can write 

p(t) = l -(I + p(t)-t) , (23) 
H = \h-B, 

and the evolution of the neutrino ensemble is determined by a precession-like equation of 
the three-vector p(t) 

f t =p(t)xh. (24) 

In this formalism attenuation effects due to CC and NC interactions can be introduced by 
relaxing the condition Tr(p) = 1. In this case 

p(f) = ^(Pb(*) , (25) 



and 

dp{E,t) 
dt 

where we have explicitly exhibited the energy dependence and 



= -i[K(E),p(E, t)] - 2X a \ E t) {n«, P(E, t)} , (26) 

--v lilt V ? / 



[Xl t (E, t)]- 1 = [\q C (E, t)]- 1 + [X NC (E, t)]- 1 , 

[Ag c (E,t)]- 1 = n T (a;)aS c ( J E), (27) 
[ANc^t)]- 1 = n T (x) a NC (E) . (28) 

tit(x) is the number density of nucleons at the point x = ct. <Jq C (E) is the cross section for 
CC interaction, v a + iV — > l a +X, and <7^c(E) is the cross section for v a + — > v a + X 
which is flavour independent. Thus we obtain four equations that describe the evolution 
of the neutrino system because one has to take into account both the flavour precession of 
the vector p(E, t) as well as the neutrino intensity attenuation encrypted in the evolution of 
p (E,t). 

v T regeneration and neutrino energy degradation can be accounted for by coupling these 
equations to the shower equations for the r flux, F T (E T ,t) (we denote by F the differ- 
ential fluxes d(f) / (dE d cos 9)) . For convenience we define the neutrino flux density matrix 
F v {E,x) = F Vfi (E ', xq) p(E , x = ct) where F u ^(E,x ) is the initial neutrino flux. The equa- 
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tions can be written as: 



FUKyX /JMEM dE , 



oo 



Anc(^,x) uy ' dE u 



dF T (E T ,t) 



-F T (E T ,x) 



dx Kec( E r,X 

-II^ ) ^ F ^ d - Es r dE - (30) 

^ ec {E T , x) = r ) T cT T . r T is the r lifetime and 7 T = E T /m T is its gamma factor. 
We have calculated the CC and NC distributions 

dNvc{E' u ,E u ) 1 da NC {E' u ,E u ) 



dE v a NC (E' v ) dE u 

dN cc (E„,E T ) _ 1 da T cc (E u ,E T ) 



dE T aQ C (E v ) dE T 



(31) 



using the MRST-g parton distributions |42j, and we have taken the r decay distribution 



dNdc j^ Ev) from Ref. |3J and djVde g;'^ } from Ref. 

The third term in Eq. represents the neutrino regeneration by NC interactions and 
the fourth term represents the contribution from v T regeneration, v T — > r~ — > v T , describing 
the energy degradation in the process. The secondary flux from v T regeneration, v T — > 
r + — ► z/ r fi + Vp, is described by the last term where we denote by over-bar the energies and 
fluxes of the r + . Br M = 0.18 is the branching ratio for this decay. In Eq. ()30|) the first term 
gives the loss of taus due to decay and the last term gives the r generation due to CC v T 
interactions. In writing these equations we have neglected the tau energy loss, which is only 
relevant at much higher energies. 

An equivalent set of equations can be written for the antineutrino flux density matrix 
and the for the r + flux. Both sets of equations are coupled due to the secondary neutrino 
flux term. 

We solve this set of ten coupled evolution equations that describe propagation through 
the Earth numerically using the matter density profile of the Preliminary Reference Earth 
Model and obtain the neutrino fluxes in the vicinity of the detector from 

d f; a f' °l = Tc[F v (E, L = 2i?cos 9) U a ] . (32) 
aE a cos a 

In Fig. El we illustrate the interplay between the different terms in Eqs. (129)1 and (|30jl. 
The figure covers the example of VLI-induced oscillations with 5c/ c = 10~ 27 and max- 
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imal £ v u mixing. The upper panels show the final and v T fluxes for vertically up- 
going neutrinos after traveling the full length of the Earth for the initial conditions 
d^(v^)o/dE v = d§(v^)Q/dE u oc E^ 1 and d<&{y T ) /dE v = d$(i/ T ) /dE u = 0. 

The figure illustrates that the attenuation in the Earth suppresses the neutrino fluxes at 
higher energies. The effect of the attenuation in the absence of oscillations is given by the 
dotted thin line in the left panel. Even in the presence of ocillations this effect can be well 
described by an overall exponential suppression [3^, |43| both for i/^'s and the oscillated 
i/ T 's. In other words, we closely reproduce the curve for "oscillation + attenuation" simply 
by multiplying the initial flux by the oscillation probability and an exponential damping 
factor: 

#„„(£, 0, L = 2Rcos6) d<j> Vllfi (E , 9) a 

dEd^d = dEd cose P ^ L = 2RcoS ^ eM-X(e)(a NC (E)+a cc (E))] . 

(33) 

where X(9) is the column density of the Earth. 

The main effect of energy degradation by NC interactions (the third term in Eq. (|29j) ) 
that is not accounted for in the approximation of Eq.(|33|) is the increase of the flux in 
the oscillation minima (the flux does not vanish in the minimum) because higher energy 
neutrinos end up with lower energy as a consequence of the NC interactions. The difference 
between the dash-dotted line and the dashed line is due to the interplay between the v T 
regeneration effect (fourth term in Eq. (|29|)) and the flavour oscillations. As a consequence 
of the first effect, we see in the right upper panel that the v T flux is enhanced because of 
the regeneration of higher energy i/ r 's, v T {E) — > t~ — > v T (E' < E), that originated from the 
oscillation of higher energies z/^'s. In turn this excess of v T 's produces an excess of z^'s after 
oscillation which is seen as the difference between the dashed curve and the dash-dotted 
curve in the left upper panel. Finally the secondary effect of v T regeneration (last term in 
Eq. (|2l?jl ). v T {E) — > r + — > n + v T v^iE' < E), results into the larger flux (seen in the left 
upper panel as the difference between the dashed and the thick full lines). This, in turn, 
gives an enhancement in the v T flux after oscillations as seen in the right upper panel. 

The lower panels show the final and v T fluxes for an atmospheric-like energy spectrum 
d^(i/fj,) /dE v = d^(Pfj)o/dE v oc E~ 3 and d$(is T ) /dE u = d$(is T ) /dE u = 0. In this case 
all regeneration effects are suppressed. Regeneration effects result in the degradation of the 
neutrino energy and the more steeply falling the neutrino energy spectrum, the smaller the 
contribution to the total flux. Therefore, in this case, the final fluxes can be relatively well 
described by the approximation in Eq. (j33"|) . 

IV. EXAMPLE OF PHYSICS REACH: VLI-INDUCED OSCILLATIONS 

Neutrino oscillations introduced by NP effects result in an energy dependent distortion of 
the zenith angle distribution of atmospheric muon events. We quantify this effect in IceCube 
by evaluating the expected angular and E^ n distributions in the detector using Eq. (J2J in 
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FIG. 4: Vertically upgoing neutrinos after traveling the full length of the Earth taking into account 
the effects due to VLI oscillations, attenuation in the Earth, v T regeneration and secondary v T 
regeneration (see text for details). 

conjunction with (and v^) fluxes obtained after evolution in the Earth for different sets 
of NP oscillation parameters. 

Together with ^-induced muon events, oscillations also generate n events from the CC 
interactions of the v T flux which reaches the detector producing a r that subsequently decays 
as t — > [w^Vt and produces a \i in the detector. Using the techniques discussed in the 
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previous sections we compute the number of z/ r -induced muon events as 

/l poo poo poo poo poo 

d cos 9 I dl dEl n / dEl I dE T / dE v (34) 
1 Jlrnin Jm u J E^ u J E° J E T 

d ^ -(E v , cos 9)^(E l/ , E T ) n T ^(E T , E°)F(E% I) A° e 



dE u dcos 9 ^ ^ v, ~dR^^ 101 ~dE^ ^ ' 

where dNdcc j^' E ^ can b e found in Ref. j^J. Equivalently we compute the number of z/ T - 
induced muon events. 

For illustration we concentrate on oscillations resulting from VLI that lead to an oscilla- 
tion wavelength inversely proportional to the neutrino energy. The results can be directly 
applied to oscillations due to VEP. 

We show in Fig. El the zenith angle distributions for muon induced events for different 
values of the VLI parameter 6c/ c and maximal mixing £ v u = 7r/4 for different threshold 
energy E^ n > -^threshold normalized to the expectations for pure Am 2 oscillations. The 
full lines include both the ^-induced events (Eq.(J2J)) and z^-induced events (Eq.(j2SJ)) while 
the last ones are not included in the dashed curves. We see that for a given value of 
Sc/c there is a range of energy for which the angular distortion is maximal. Above that 
energy, the oscillations average out and result in a constant suppression of the number of 
events. Inclusion of the ^-induced events events leads to an overall increase of the event 
rate but slightly reduces the angular distortion (see also Fig. EJ) as a consequence of the 
"anti-oscillations" of the z/ r 's as compared to the z/ M 's. 

In order to quantify the energy-dependent angular distortion we define the vertical-to- 
horizontal double ratio 

Nf{El m '\ -0.6 < cos9 < -0.2) 
M " ' Pver 11 N vh {E fm '\-l< cos9 < -0.6) V ; 



N™- vl \El m <\ -1< cos# < -0.6) 

where by E^ m,i we denote integration in an energy bin of width 0.2 log 10 (£^ m '*) using that 
IceCube measures energy to 20% in log 10 E for muons. 

In what follows we will use the double ratio in Eq. (|H5|) as the observable to determine 
the sensitivity of IceCube to NP-induced oscillations. We have chosen a double ratio to 
eliminate uncertainties associated with the overall normalization of the atmospheric fluxes 
at high energies. 

Also, in the definition of the double ratio we have conservatively included only events well 
below the horizon cos# < —0.2 to avoid the possible contamination from missreconstructed 
atmospheric muons which can still survive after level 2 cuts in the angular bins closer to the 



horizon 40 ] . 



In Fig.lHlwe plot the expected value of this ratio for different values of 5c/ ' c. As mentioned 
above, IceCube measures energy to 20% in log 10 E for muons. Accordingly, we have divided 
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FIG. 5: Zenith angle distributions for muon induced events for different values of the VLI parame- 
ter Sc/c and maximal mixing £ v n = ir/4 for different threshold energy E^ n > ^threshold normalized 
to the expectations for pure Am 2 oscillations . The dashed line includes only the f^-induced muon 
events and the full line includes both the z^-induced and v T - induced muon events. 
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FIG. 6: The predicted horizontal-to- vertical double ratio in Ea. ()35|) for different values of 5c/ c. 
The data points in the figure show the expected statistical error corresponding to the observation 
of no NP effects in 10 years of IceCube. 



the data in 16 E^ m bins: 15 bins between 10 2 and 10 5 GeV and one containing all events 
above 10 5 GeV. In the figure the full lines include both the ^-induced events (Eq.(J2J)) and 
z/ r -induced events (Eq.()35p) while the last ones are not included in the dashed curves. As 
described above, the net result of including the z/ r -induced events is a slight decrease of 
the maximum expected value of the double ratio. The data points in the figure show the 
expected statistical error corresponding to the observation of no NP effects in 10 years of 
IceCube. 

In order to estimate the expected sensitivity we assume that no NP effect is observed and 
define a simple x 2 function as 



^ (R h/v (E f m >\5c/c,Z vli )-l) 2 
X {Sc/c, £ v ii) = 2^ 

i=i 



a: 



(36) 



where <7 s t a t,i is computed from the expected number of events in the absence of NP effects 
(see Table |TJ). 

We show in Fig. the sensitivity limits in the [Sc/c, £ v h] -plane at 90, 95, 99 and 3 o 
CL obtained from the condition x 2 (<5c/c, £ v ii) < Xmax (CL, 2dof). In order to estimate the 
uncertainty associated with the poorly known prompt neutrino fluxes we show in the figure 
the results obtained using the RQPM model (filled regions) and the TIG model (full lines). 
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TABLE I: Number of expected atmospheric z^-induced muon events in 10 years of IceCube oper- 
aton in the different energy bins and angular bins used in the analysis, assuming no NP effect is 
observed. 



The difference is about 50% in the strongest bound on 5c/ ' c. 

The figure illustrates the improvement on the present bounds by more than two orders of 
magnitude even within the context of this very conservative analysis. The loss of sensitivity 
at large 5c/ c is a consequence of the use of a double ratio as an observable. Such an 
observable is insensitive to NP effects if 5c/ c is large enough for the oscillations to be always 
averaged leading only to an overall suppression. 

When data becomes available a more realistic analysis is likely to lead to a further im- 
provement of the sensitivity. 

V. SUMMARY 

In this paper we have investigated the physics that can be probed with the high-statistics 
high-energy data on atmospheric neutrinos which will be collected by the IceCube detec- 
tor. In order to do so first we have developed a semianalytical simulation of the detector 
performance. In particular, we present in Eqs. (J2J)- (jlU)) the parametrization of the effective 
area of the detector which correctly reproduces the results of the detailed experimental MC 
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FIG. 7: Sensitivity limits in the 5c/c,£ v n at 90, 95, 99 and 3 a CL. The hatched area in the upper 
right corner is the present 3cr bound from the analysis of SK data in Ref. 



for the response of the IceCube detector after events that are not neutrinos have been re- 
jected using the quality cuts referred to as level 2 cuts. We conclude that in 10 years of 
operation IceCube will collect more than 700 thousand atmospheric neutrino events with 
energies Ejj 111 > 100 GeV which offer a unique opportunity to test new physics mechanisms 
for leptonic flavour mixing which are not suppressed at high energy. In general these effects 
are expected to induce an energy dependent angular distortion of the events. 

Next, because of the relatively high energy of the neutrino sample, NP induced flavour 
oscillations, propagation in the Earth, regeneration of neutrinos due to r decay must be 
treated in a consistent way. In Sec. IHII we have presented the corresponding evolution equa- 
tions. We conclude that for steeply falling neutrino energy spectra, such as the atmospheric 
neutrino one, the dominant effect together with flavour oscillations is the attenuation of the 
oscillation amplitude due to inelastic CC and NC interactions of the neutrinos in the Earth in 
conjunction with the production z/ r -induced muon events due to the chain v T — > r — > ^jlv^Vt 
in the vicinity of the detector. iv-induced muon events can increase the event sample by at 
most O(10%). 

Finally we have applied these results to realistically evaluate the reach of IceCube in 
studying physics beyond conventional neutrino oscillations induced by violation of Lorentz 
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invariance and/or the equivalence principle. In Fig. [7|we show how even with a very con- 
servative analysis the range of testable sizes of these effects can be easily extended by more 
than two orders of magnitude. The methods developed are readily applicable to probe 
speculations on other non-conventional physics associated with neutrinos. 
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